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1. INTRODUCTION 

A greenhouse is a structure that can be perfectly enclosed to protect plants and to promote the 
growth of crops (vegetables, flowers, etc.) in creating more favorable climatic conditions than the local 
climate. It aims to remove food crops from climatic elements for better management of plant needs and to 
accelerate growth or produce them regardless of the seasons. 

The main asset of greenhouses is crop protection. Indeed, the greenhouse effectively protects against 
changes in climatic conditions, such as cold or hail. By monitoring the temperature, relative humidity, CO2 
level. And good aeriation with a combination of a roof vent, front doors and fans [1]. 


1.1 Background 

The greenhouse is a multivariable non-linear system. Thus, with the current control methods used, 
for the synthesis of the parameters of the dynamic model of the greenhouse, mathematical models likely to 
represent the greenhouse will be necessary to identify [2]. 


1.2 The Problem 

The problem of optimal control when applied to coupled multivariable systems is the adaptability 
between its control strategy and the internal model of the controller process, which is generally very time- 
evolving [3-4]. 


1.3 Proposed Solution 
The linear Gaussian linear method known as LQG will be integrated with an adaptive strategy 
control using an on-line identification approach. The LQG control is a technique that is widely used and 
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recently implemented in various industrial applications [5]. It makes it possible to calculate the gain of a 
control by return of state in a particular concern to reduce the white noises. However, the LQG control 
combines an LQ controller (Linear Quadratic) and a Kalman estimator that can be calculated independently 
according to the separation principle. 

The paper is organized as follows: The first part is devoted to the modeling of the greenhouse and its 
dynamic behavior using energetic transfers to its indoor climate. In the second part, we present the linear 
quadratic Gaussian (LQG) optimal control technique. The adaptive control scheme combined with parametric 
identification of the process is developed in the third part. In the fourth part, we present the results of 
simulations obtained by the application of the LQG controller on the climate control of the greenhouse. A 
general conclusion will summarize the work done, as well as the prospects for possible improvements. 


2. RESEARCH METHOD 

The models based on energy storage elements give qualitatively good results when idealized input 
signals are applied. The greenhous [6]. is a multi-input and multi-output (MIMO) system which is equipped 
with several sensors and actuators as show in Figure 1. 


To Ho Sr Sw 


Rc 
Greenhouse 


Figure 1. Schematic Diagram of Controlled Greenhouse 


In a first approximation the microclimate greenhouse shown on Figure | is described by: 
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with 


Rc: heating energy applied to the plant (Kw). 
Vt: ventilation angle outside the greenhouse (°C). 
To;Ho: air temperature and relative humidity outside the greenhouse (°C, %100). 
Sr: Solar radiation (W / m). 
Ti;Hi: air temperature and relative humidity inside the greenhouse (°C, %100). 
Sw: wind speed outside the greenhouse (m/s). 
We consider that the coefficients a;;,bjjand dj; are slowly varying in time,then the equations 
overhead can be written as: 


ce +1) = AX(k) + BU(k) + DNd(k) 


¥(k) = CX(k) a) 


3. LINEAR QUADRATIC GAUSSIAN (LQG) CONTROL 
3.1 Design of the Controller 

The objective of LQG controller seen in Figure 2 is to generate a suite of commands u which allow 
the output to follow a desired behavior of the greenhouse determined by the needs of the plant [7-8]. These 
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commands are calculated according to a criterion optimized, from the state space representing the 
microclimate greenhouse. Consider the following stochastic system: 








Figure 2. Structure of the LQG Control Method. 


x{t) = Ax(t) + Bu(t) + Le (t) ’ 
5 (t) = Cx (t) + Du(t) + Ke (t) (4) 


where u is the input vector of the process, y is the output, L and K represents Gaussian sequences of zero 
mean with covariance, and x (k) is the state vector. The objective is to predict, from the initial state, a 
sequence of inputs u that allows the output y of the process to follow a reference trajectory set by the user. 
The quadratic cost criterion which we seek to minimise in the LQ problem is given by 


J= So x™® Qx(t) + u™@Ru(t)]dt (5) 


Q and R are non-negative definite symmetric matrices, while the state observer is a filter to rebuild 
the system status from its dynamic model, is usually represented in the form of state space. The output of this 
filter is a close vector of the system state vector. Consider the following deterministic system (without 
random noise) [9]. 


f (t) = Ax(t) + Bu(t) (6) 
y(t) = Cx (t) 


From which we built the observer following status: 


X(t) = AX (t) + Bu(t)4+ Lv (t) 
H(t) = CX) ee) 
v(t) =y@) -I@ 


The error between the state vectors is written as follows: 


(FO = x(t) — £(t) 


é() = 4% () -¥ (1) (8) 
By replacing the state vectors by their expression of the Equation 7, we get this 

é6(t) = Ax (t) + Bu (t) — [AZ (t) + Bu (t) + L(y() -9@)] (9) 
Then we have 

é(t) = [A — LC] e (t) (10) 
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If the eigenvalues of matrix [A — LC] are all in the complex left half-plane, the system is 
asymptotically stable and the error vector between the state and x the state X tends to zero exponentially [10]. 
Thus, we can use a state observer and place an order to return state using the observed state such as 


x(t) = Ax (t) — BKR(t) + w(t) 
le (t) = AR(t) — [BK + LC]£(t) + Lv(t) + LCx (t) ey 
Which can be written as: 

x(t) = Ax (t) — BK [x (t)- e@]+ w(t) (12) 
é(t) = [A- LC]e(t) + w(t)- Lv(t) (13) 
Or, gathered in the following matrix form: 

Eels oe acidleeltl Ste (1) 
The quadratic cost criterion objective [11] can be rewritten as 

Jugc = ECSy [x"(t) Qx(t) + u™ (t)Ru(t)]dt) (15) 
Then, 

Joa = ES) (& + 8)" (t) QR + eC) + u?(t)Ru(t)] at) (16) 
We get, 

Jrac = Sy (8) Qx(t) + u(t)Ru(t)]dt + 2 f> [7 (0) QE(e)]dt + E Jy [é" (t) Qelat (17) 
e is acentred random variable, then E(e) = 0, it comes, 

Iroc =Jug t ES, [67 (t) Qelat (18) 


With, Ji9 is a type of criterion LQ on the observer state vector & of the plant. the observer must be 


designed so that the quantity 1h [€7(t) Qe] is minimal, this observer is called optimal Kalman [12]. 


3.2. Formulation of the Discrete LQG Controller 
The system is described by the following discrete equation: 


{20 = Ax (k) + Bu(k) + Le (k) (19) 


y (kK) = Cx (kK) + Du(k) + Ke (k) 


The structure of the control u is conditioned by using the criterion appeared in the equation 16. 


J = limy-soo 7 Dho((K) — y-(K))” + Auch)? 
(20) 


With N, the horizon of LQG control yr is as set point desired signal , and / is the weighting factor 
for modulating the u input sequences [13]. Hence the optimal control minimise this criterion is 


u(k) = —Lx(k) + r(k) (21) 
With, 
L(k) = ((A + BTR (k)B)"1B'V (k)) (22) 
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V (k) = (A — BL(k))" — Cy,(k) (23) 
R(k) = ATR(k —1)A—A™R(k — 1)B(A + BTR(K — 1)A+C"C) (24) 


With R(O) > 0 . The solution of algebraic Riccati equation, its existence depends to two conditions: 
- The pair (A, B) is to be stabilisable. This condition ensures the stability of the LQG control. 
- The pair (C, A) should be observable. This condition ensures the stability of the estimator. 


x(k +1) = Ax(k) + Bu(k) + K (k) (vy (ke) — Cx(k)) (25) 


The current state estimate x (kK) directly computed the equation above, With K (k) is the gain of the 
Kalman filter adaptive algorithm. The stationary model previously used do not take into account the 
nonlinearities. This model often underestimates the order process. In addition, the structure of the process 
varies slowly over time (growth of the vegetation, walls aging). Therefore, we propose to use an adaptive 
model. 


3.3. Parametric Identification Method 

The adaptive form of the control algorithm [14-15]. described above consists of using the estimate 
of the plant model parameters instead of the unidentified real values in the control law. To estimate the 
parameters of the model described above, the recursive least square (FFRLS) algorithm with forgotten factor 
has been used [16]. This method uses the model given in equation 


Y(k) =O" (K)e(k) (26) 
Where 6 is the vector of the unknown parameters defined as 
67(k) = [a11(k), ...,d2q(k)] (27) 


In Equation 26, @(k) is a regression vector partly consisting of measured input/output variables and 
is defined as 


g(k) = [Ti Hi Rc Vt To Ho Sr Sw | (28) 
In case the parameters of the model change slowly, we use the method of recursive least squares 


with exponential forgetting factor. and then try to calculate the parameters that minimize a weighted criterion 
by a power of a coefficient A [17]: 


6(k) = Arg min YL, aK (yk) — 97 (k).6(kK))? (29) 


Thus the oldest residues are less important. And by a short demonstration, the following 
identification algorithm is given [18-19]: 


8(k) = O(k — 1) + K(k)[y(k) — 97 (k).8(k — 1] (30) 
K(k) = P(k— De + oT (K)P(k- Dek)? (31) 
P(k) = = (P(k — 1) — K(k)g" (k)) (32) 


For, the gain matrix P(t) will not converge to zero as time goes to infinity, hence the algorithm will 
be able to update the parameters when the system is time varying. In practice, is a design parameter a low 
value of gives a fast tracking of time varying parameters but a high noise sensitivity [20]. Otherwise, one will 
choose A(k) < 1 (typically between 0.98 and 0.995. The value of A depends on the dynamics of changing 
parameters of the system. 


4. SIMULATION RESULTS 
The objective of this study is to show how adaptive LQG implementation can cope with 
Multivariable systems. The LQG design is applied to the microclimate greenhouse plant with changes in 
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external disturbances. The Simulation of the plant with the proposed adaptive control scheme are 
implemented in real time process simulation with the measured weather values of the air temperature and 
relative humidity and other disturbances (solar radiation, wind speed, external temperature and relative 
humidity) a Gaussian noise with constant variance is added to disturbances inputs of a greenhouse. From the 
measured input output signals, the identified model of the microclimate greenhouse that give the best 
behavior discrete model is: 


0.4752 0.0878 pa Pi 0.0165 
0.7879 0.5004 0.1084 0.0704 


_ 0.3450 —0.0385 0.0138 0.2012 
—0.7095 0.3639 —0.0081 0.2765 


Where the sample time is T = 5 seconds, the reference to be track is chosen as a square wave. After 
several checks, the recursive identification is made via the FFRLS algorithm with the flowing consideration 
as stated in [21] 

- forgetting factors are 1 = 0.97 

- the initial adaptation gain of P(O) = 10*I 

- the initial parameter value is 67 =[0.01, 0.01, 0.01]. 

- the following weighting matrices for the Kalman filter Q = [100 0; 0 50]; R = [10 0; 0 15], 

The implementation of the adaptive FFRLS algorithm and the LQG multivariable control for 
monitoring the temperature and humidity inside the greenhouse, using the MATLAB programming code, 
gave us the following results. The outputs predicted by the adaptive identification method rigorously 
converge with the outputs of the model simulating the thermal behavior of the greenhouse, as can be seen 
from Figure 3. 


Evolution of real and estimated temperature 
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Figure 3. Measured and Simulated Responses of the Greenhouse Air Temperature and Relative Humidity 


Time-based Evolution of the Adjustable Parameters of the microclimate green house shown in 
Figure 4. The Figure 4 has allowed us to highlight that: 

Firstly, the recursive FFRLS algorithm for estimating the observation vector based on the square 
error update produces more reliable estimates of the greenhouse parameters in the presence of the other 
disturbances. 

Secondly, this method using the recursive technique with a forgetting factor has a fast convergence 
dynamic that makes it easily applicable on changing systems. The recursive least squares with forgotten 
factor method presented in this paper gives unbiased estimates only for ARMAX models. On the other hand, 
this method, does not use any information a priori on the measurement noise, and we have tested that if the 
noise is not zero average value. the estimation of the parameters is biased. 
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Figure 4. Time-based Evolution of the Adjustable Parameters of the Microclimate Green House (see online 
version for colours) 


To illustrate the applicability and efficiency of LQG control technique on the greenhouse process. 
Figure 5 show the real-time behavior of temperature T(°C) and humidity H(°C) respectively. These results 
once again confirm the robustness of the LQG control to be able to cancel almost completely the effect of the 
strong disturbances added to the thermal process of the greenhouse compared to the standard monovariable 
control. 
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Figure 5. Responses of Greenhouse Air Temperatures and Relative Humidity using the Adaptive LQG 
Controller 
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With the observation of the Figure 6 , we note that in spite of the presence of the coupling between 
the variables of the greenhouse, the change of setpoint carried out in real time on both the temperature (high) 
and the second output which is the humidity of the air (bottom) has less impact on both outputs. The 
commands applied to the system are given in Figure 6. A s in simulation mode, the influence of the 
multivariable coupling is clearly felt on the first heating control at the time of the change of the set point 
applied to the second air humidity output H (%). The first control reacts so that the temperature of the air T (° 
C) is not disturbed by the change in behavior of the greenhouse that affects humidity H (%). 





Input Control Heating & Ventilation 














Figure 6. The Input Control of the Greenhouse with the Adaptive Multivariable LQG Controller 


5. CONCLUSION 

In this paper, a multivariable adaptive LQG controller combined with the FFRLS parameter 
estimation method was applied to driving the microclimate greenhouse air temperature and relative humidity 
with external disturbances. The simulation results show that the proposed scheme is able to sufficiently 
control the plant, and conserves a high level of performances, in terms of pursuing, response time and 
disturbance rejection. 
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